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I.  INTRODUCTION 


The  objective  of  this  investigation  is  the  development  of  a  theoretical  model 
for  propagation  of  Extremely  High  Frequency  (F.HF)  radio  waves  in  a  strongly  scat¬ 
tering  medium  such  as  a  forest.  The  forest  is  described  as  a  random  medium  which, 
in  first  approximation,  may  be  regarded  as  statistically  homogeneous;  and  wave  prop¬ 
agation  in  this  medium  is  characterized  by  the  theory  of  radiative  energy  transfer 
(transport  theory).  This  theory  permits  one  to  take  multiple  scattering  effects 
fully  into  account  while  interference  effects  are  neglected.  Experimental  evidence 
confirms  that  the  coherence  length  of  mm  waves  in  vegetation  is  short  so  that  the 
suppression  of  interference  effects  is  an  acceptable  approximation. 

All  scatter  objects  in  a  forest  environment,  including  tree  trunks,  branches, 
and  leaves,  or  pine  needles,  have  large  dimensions  compared  to  a  wavelength  in  the 
nun-wave  region.  As  a  consequence,  the  forest  medium  will  produce  strong  forward 
scattering  and  its  scatter  characteristic  (phase  function)  may  be  represented  by  a 
pronounced  forward  lobe  with  an  isotropic  background.  Assuming  this  type  of  phase 
function,  the  medium  is  characterized  by  four  parameters:  the  absorption  cross 
section  per  volume  o^,  the  (total)  scatter  cross  section  per  volume  o^,  the  beam 
width  Ays  of  the  forward  scatter  lobe,  and  the  ratio  a  of  the  forward  scattered 
power  and  the  total  scattered  power. 

These  four  quantities  are  determined  by  the  macroscopic  and  microscopic 
structure  parameters  of  the  forest  and,  in  principle,  can  be  derived  from  these 
parameters.  No  attempt  is  made  here  to  establish  this  relation  theoretically. 
Considering  the  structural  complexity  of  the  forest  environment,  this  would  appear 
to  be  a  hopeless  task.  Rather,  the  transport  equation  is  formulated  and  solved  in 
general  terms  and  the  dependency  of  the  solutions  on  the  parameters  o^,  o^,  Ay^  and 
a  is  established.  Comparison  with  experimental  results  then  permits  one  to  deter¬ 
mine  estimates  for  these  parameters  by  adjusting  them  such  that  the  theoretical 
curves  agree  with  the  experimental  data.  At  the  same  time,  the  theory  serves  to 
explain  and  interpret  the  experimental  results.  Theory  and  experiments  are  shown 
to  be  in  good  qualitative  agreement.  The  experimental  data  utilized  in  this  report 
was  obtained  by  Violette  and  Espeland  1.1,  21.  Achieving  good  quantitative  agree¬ 
ment,  however,  is  not  without  difficulties.  A  larger  measured  data  base  and  a  re¬ 
finement  of  the  theory  will  probably  be  required. 


Propagation  characteristics  to  bo  determined  include  the  range  dependence  of 
both  the  coherent  and  incoherent  field  components  of  mm-wave  beams  transmitted 
through  woods  and  forests,  and  the  directional  spectrum  of  the  incoherent  component 
(beam  broadening)  as  a  function  of  vegetation  depth  and  density.  The  theory  shows, 
in  particular,  that  the  range  dependence  is  not  determined  by  a  simple  exponential 
decrease  in  signal  strength  at  constant  attenuation  rate  but  a  high  attenuation 
rate  at  short  distances  into  woods  gradually  transforms  into  a  much  lower  attenu¬ 
ation  rate  at  large  distances.  This  complex  range  dependence  has  been  confirmed 
by  the  experiments  and  is  explained  by  the  interplay  of  the  coherent  (direct  path) 
component,  which  dominates  at  short  distances  into  woods  but  is  strongly  attenu¬ 
ated,  and  the  incoherent  (multiply-scattered)  component,  which  is  less  attenuated 
and  takes  over  at  large  distances.  The  transition  is  accompanied  by  substantial 
beam  broadening.  The  comparison  of  theory  and  experiments  indicates,  furthermore, 
that  the  parameter  a  is  close  to  unity  and  that  o<-  >>  In  other  words,  forward 
scattering  seems  to  be  strong  and  scattering  seems  to  dominate  over  absorption. 


J .  Fundamental  $ 

The  electromagnetic  field  in  a  random  medium  of  scatterers  can  be  split  into 
two  parts,  the  coherent  component  and  the  incoherent  component.  The  coherent,  com¬ 
ponent  is  the  average  field  which  at  any  point,  in  space  has  a  well-defined  direction 
of  propagation  and  a  well-defined  polarization.  This  component  decreases  due  to 
both  absorption  and  scatter,  and  its  attenuation  rate  is  comparatively  high.  The 
incoherent  component  is  the  zero-mean  field.  It  is  generated  by  the  scattering  of 
the  coherent  component  and  consists  of  many  wave  trains  propagating  in  various  di¬ 
rections.  Hence,  at  any  point  in  space,  the  incoherent  component  does  not  have  a 
unique  direction  of  propagation  but  is  characterized  by  a  directional  spectrum.  In 
transport  theory,  it  is  assumed  that  the  constituent  wave  trains  are  uncorrelated 
in  phase  and  add  in  power.  The  incoherent  component  decreases  at  a  slower  rate 
than  the  coherent  component  since  it  is  attenuated  due  to  absorption  only  (though 
over  an  increased  path  length)  while  it  scatters  into  itself.  Scattering  of  the 
incoherent  component  reproduces  the  incoherent  component..  In  general,  this  field 
component  is  partially  depolarized. 

While  the  coherent  component  will  dominate  at  short  distances  from  a  coherent 
source,  the  incoherent  component  takes  over  at  large  distances,  leading  to  a  re¬ 
duced  attenuation  rate,  significant  beam  broadening,  and  depolarization.  The  co¬ 
herent  and  incoherent,  field  components  are  powerwise  orthogonal.  If  E  is  the  elec¬ 
tric  field  strength  in  a  random  medium,  then  the  average  field  intensity  is  given 
by 

-  <ir><r* >  +  <aF’Af.*>  _  _ 

with  AE  -  E  -  <E> 


In  other  words,  the  total  intensity  is  the  sum  of  the  intensities  of  the  coherent 
and  incoherent  field  components. 

The  fundamental  quantity  used  in  transport,  theory  is  the  specific  intensity 
I(r,s)  which  is  defined  as  the  power  per  unit  area  and  per  unit  solid  angle  propa¬ 
gating  at  point  F  into  the  direction  s: 

■  IF# 


in 


Here,  da  is  the  area  of  the  surface  element  (at  point  F)  normal  to  the  unit 
vector  s,  and  di2  is  the  differential  of  the  solid  angle  ST  centered  about  f;  see 


Fig.  1.  Definition  (1)  applies  to  time  harmonic  fieUls  in  a  time  independent  ran¬ 
dom  medium,  the  case  considered  hero.  (Wind  induced  motion  in  woods  and  forests  is 
negligibly  slow  compared  to  the  travel  times  of  EM  signals.)  The  specific  inten¬ 
sity  can  be  used  to  characterize  both  the  coherent  and  incoherent  field  components. 
In  the  case  of  the  coherent  component,  which  has  a  unique  direction  of  propagation, 
I  takes  the  form  of  a  6- function  in  s'. 

The  transport  equation  [3],  the  basic  equation  of  the  theory  of  radiative 
energy  transfer,  is  obtained  by  formulating  the  power  conservation  theorem  in  terms 
of  the  specific  intensity  I.  Consider  the  cylindrical  volume  dV  =  da-dl  of  Fig.  2, 
which  is  centered  at  point  r  and  whose  axis  has  the  direction  of  the  unit  vector  s. 
Radi  a  t i on  with  the  specific  intensity  I  (r,s)  passing  through  this  volume  will  vary 
because  of  two  effects: 

(1)  Absorption  and  scattering  by  the  random  medium  will  reduce  the 
power  propagating  in  the  direction  s'.  The  specific  intensity 
will  decrease  by  the  amount 

dl>  -  -(uA  i  os)  Kr.s)  dl  (2) 


(2) 


where  is  the  absorption  cross  section  per  unit  volume  of 
the  random  niodiuin  and  is  the  (total  )  scatter  cross  section 
per  volume.  In  general,  both  o^  and  depend  on  r. 


A  fraction  of  the  power  propagating  in  other  directions  s' 
will  be  redirected  into  the  direction  s  by  scattering  with-, 
in  the  volume  element  dV.  This  effect  will  increase  I  (Y\T) 
by  an  amount  which  is  determined  by  the  scatter  characteris¬ 
tic*  of  the  medium  p  (?,?'')•  The  increase  of  I(F,s)  is  given 
by 

d l?  -  a<.  IJ p(s.s')  I(r,s'D  da'  dl 


where  the  integration  is  performed  over  all  spatial 
directions  z'  and  p  is  assumed  normalized  such  that 
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(4) 


*Ir.  transport  theory  p(s.s')  is  usually  called  the  "phase  function"  of  the  scatter 
medium,  though  it  is  not  related  in  any  way  to  the  phases  of  the  electric  or  mag¬ 
netic  field  strengths. 


d 


Similar  to  and  o<,,  the  scatter  characteristic,  in  general, 
will  be  a  function  of  r.  Reciprocity  requires  that  p  is  sym¬ 


metric  in  s  and  s  : 


p(s.s')  =  p ( s  , s ) 


The  total  variation  of  the  specific  intensity  I  over  the  length  dl  of  the 
volume  element  dV  is  the  sum  of  dl^  and  dl^-. 

dl  ^  [s-grad  I(F,T)]  dl  =  d I ^  +  dl? 

Hence,  with  Equations  (2)  and  (3): 


•grad  l(r,s)  +  (a.  +  o<.)  I(r,s) 


Jf  "<*• 


s')  I(r,$')  dft' 


Equation  (6)  is  the  transport  equation  in  its  scalar  formulation.  A  forcing 
term  which  usually  appears  on  the  right  side  of  this  equation  is  omitted  here  since 
we  will  assume  that  all  sources  are  located  outside  the  random  medium.  More  pre¬ 
cisely,  we  shall  assume  that  the  transmit  antenna  is  situated  outside  the  forest 
at  a  certain  distance  from  the  forest-to-air  interface  while  the  receive  antenna  is 
placed  within  the  forested  region.  Reciprocity  insures  that  the  same  results  will 
be  obtained  when  transmitter  and  receiver  are  interchanged. 

2.  Formulation  of  the  Problem 

In  analyzing  nun-wave  propagation  in  woods  and  forests,  we  use  a  simple  model 
which  is  illustrated  in  Fig.  3.  A  forest  halfspace  is  separated  by  a  planar  in¬ 
terface  from  an  air  halfspace.  The  forest  is  characterized  as  a  statistically 
homogeneous  medium  of  random  scatterers.  In  other  words,  it  is  assumed  that 
u<.  and  p  do  not  vary  with  r  and  are  the  same  at  every  point  in  the  medium.  Further¬ 
more,  we  assume  that  the  phase  function  p  depends  only  ori  the  scatter  angle  y  - 
arc  cos  (§~-s^)  subtended  by  the  directions  s  and  s'.  Since,  in  a  forest,  all  scat¬ 
ter  elements  are  large  compared  to  nun-wavelengths,  strong  forward  scattering  will 
occur  and  d  is  assumed  to  consist  of  a  narrow  Gaussian  forward  lobe  superimposed 
over  an  isotropic  background,  i.e.. 


P(y)  -  a  (J^)2  e 


_(JL  )? 

5  +  (1-x) 


\ 


FIGURF  3:  Forest  Halfspace  (Z  >  0) 
Illuminated  by  Plane  Wave 


The  function  p(y)  is  indicated  in  Fig.  4.  The  width  of  the  forward  lobe  Ays 
is  assumed  to  be  «tt,  and  a  is  the  ratio  of  the  forward  scattered  and  the  total 
scattered  power.  It  is  easily  verified  that  p(y)  satisfies  Equation  (4). 

The  transmitter  ir.  the  air  halfspace  is  assumed  to  be  sufficiently  far  away 
from  the  edge  of  the  forest  so  that  the  field  incident  upon  the  air-to-forest  in¬ 
terface  can  oe  described  as  a  plane  wave.  This  plane  wave  may  be  impinging  nor¬ 
mally  or  at  an  oblique  angle.  The  coordinate  system  is  indicated  in  Fig.  3.  In 
particular,  Z  is  the  coordinate  normal  to  the  interface  which  is  assumed  to  coin¬ 
cide  with  the  plane  Z  =  0;  0  is  the  angle  counted  from  the  positive  Z-direction; 
and  (|>  is  the  projected  angle  in  planes  Z  -  constant. 


With  the  above  assumptions,  the  specific  intensity  I  in  the  forest  medium 
will  depend  on  the  coordinate  Z  and  the  angles  8  and  $  only,  but  it  is  independent 
of  the  cross-sectional  coordinates  x  and  y.  The  transport  equation  (6)  then  re¬ 
duces  to 


cos  8  (Z;  0 ,<J>)  +  (oA  +  os)  I(Z;  6,<f>) 

2lf  TT 

=  4“  J  J  p(Y)  HZ;  sin  tj'de'd<T 

0  o 

for  Z  >  0 

where 


(8a) 


cos  y  =  cos  (4>-<t>")  sin  8  sin  6"  +  cos  0  cos  8'  (8b) 

In  the  plane  Z  =  0,  I  is  determined  by  the  incident  plane  wave  (for  all  for- 
Honce,  we  have  the  boundary  condition 

<5(0-61 


ward  directions 


■«>•.<>.♦>  ■  sp  S«"V 


for  0  <  0  < 


ir 

7 


(9a) 


The  symbol  6  denotes  the  well  known  6-function.  S  ,  0  and  <j>  are  the 

r  r  r 

Foynting  vector  and  the  propagation  angles  of  the  incident  wave.  Finally,  since 
the  medium  is  not  only  scattering  but  also  dissipative,  I  must  vanish  at  infinity, 
i  .e. , 

I  (Z;  0,4))  •*  0  for  Z  m  (9b) 

Equations  (8)  and  (9)  provide  the  transport  theory  description  of  the  forest 
model  considered  here;  the  equations  permit  a  unique  solution  of  the  problem. 
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FIGURE  4:  Phasefunction  with  Narrow  Forward  Lobe 
and  Isotropic  Background 


Mathematically  speaking.  Equations  (8)  and  (9)  formulate  the  problem  in 
terms  of  a  homogeneous  integro-differential  equation  with  an  inhomogeneous  bound¬ 
ary  condition.  By  a  simple  transformation,  the  problem  may  be  reformulated  as  an 
inhomogeneous  equation  with  a  homogeneous  boundary  condition.  The  transformation 
is  discussed  here  since  it  provides  additional  physical  insight.  We  split  the 
specific  intensity  1  into  two  parts  [3]: 

i(z,e,4>)  =  Iri  (z;  e,t(>)  +  Id  (z;e,*)  (10) 


where  1  -  is  called  the  reduced  incident  intensity  and  I  .  the  diffuse  intensity 
ri  d 

[3]  .  is  chosen  to  satisfy  the  boundary  condition  in  the  plane  z  =  C  and  the 

transport  equation  (8a)  with  the  interaction  term  on  the  right  side  suppressed, 
i  .e. 

5 (0 “0  ) 

Iri  (0;O ,<.)  =  Sp  for  0  <  e  <  \  (Ha) 


3Iri 


cose  -T.— -  +  (ofl  +  Of.)  0 


for  z  >  0 


(lib) 


These  equations  can  be  solved  in  closed  form: 

oA  +  ot 


Vi  ■  sp  “of-fesr  «  ■  ^  S<*-V 


for  Z  >  0 


(12) 


The  diffuse  intensity  then  has  to  satisfy  the  equations; 


(  0  i  0  , 0 )  - 

0 

for  0  <  0  <  j  ,  Z  =  0 

(13a) 

Id  (z ;6 ,4>)  -*■ 

0 

for  7.  -*■  °° 

(13b) 

91  . 

37  +  (aA  + 

2ir  n 

_°s 

4tt 

//*» 

Id  (Z;  6 ',<!>')  sin  O'dG'di))' 

u  o 
0  0 


°AaS  (13c) 

a  i-Vi  2. 


for  Z  >  0 


with  cos  Yp  =  cos  sin  8  sin  8p  +  cos  e  cos  0p  (13d) 

The  second  term  (forcing  item)  on  the  right  side  of  Equation  (13c)  is  the  contri¬ 
bution  of  I  ^  to  the  interaction  integral. 

Evidently,  the  reduced  incident  intensity  Ir<  is  the  continuation  of  the 
incident  plane  wave  into  the  random  medium  where  its  magnitude  decreases  exponen¬ 
tially  with  distance  because  of  absorption  and  scattering.  Similar  to  the  inci¬ 
dent  plane  wave,  Ir^  has  a  well-defined  direction  of  propagation;  it  can  be  con¬ 
sidered  as  the  coherent  field  component,  i.e.,  as  the  specific  intensity  of  the 
average  field. 

The  diffuse  intensity  1^,  on  the  other  hand,  can  be  interpreted  as  the  in¬ 
coherent  field  component,  i.e.,  it  represents  the  specific  intensity  of  the  zero- 
mean  field.  The  incoherent  component  is  generated  by  the  scattering  of  the  co¬ 
herent  component  as  represented  by  the  inhomogeneous  term  in  Equation  (13c).  The 
energy  flow  of  1^  directed  into  the  scatter  medium  is  zero  in  the  interface  plane 
Z  =  0  where  no  scattering  of  the  coherent  component  has  taken  place  as  yet.  This 
is  formulated  by  boundary  condition  (13a). 

Before  proceeding  to  solve  Equations  (13),  the  problem  of  characterizing 
the  polarization  properties  of  electromagnetic  fields  in  random  media  is  briefly- 
discussed  here.  The  coherent  component  has  the  same  well-defined  polarization  as 
the  incident  plane  wave.  The  incoherent  component,  on  the  other  hand,  is  in  gen¬ 
eral,  partially  depolarized.  This  means  that,  in  actuality,  the  diffuse  component 
should  be  represented  in  terms  of  the  four  Stokes  parameters  and  the  phase 
function  p  in  terms  of  a  4  x  4  element  Stokes  matrix,  which  would  result  in  a  form¬ 
ulation  of  the  radiative  energy  transfer  problem  in  a  forest  in  terms  of  four 
coupled  integro-differential  equations. 

However,  if  reasonable  assumptions  are  made  concerning  the  scatter  prop¬ 
erties  of  the  forest  medium,  the  Stokes  matrix  can  be  simplified  and  the  equation 
for  the  first  Stokes  parameter  becomes  decoupled  from  the  remaining  three 
equations.  The  first  Stokes  parameter  is  of  particular  interest  since  it  repre¬ 
sents  the  sum  of  the  intensities  associated  with  two  orthogonal  polarizations, 
i.e.,  it  represents  the  total  intensity  of  the  incoherent  component  at  any  point 
in  space.  The  decoupled  equation  for  this  first  Stokes  parameter  takes  the  form 
of  the  scalar  transport  equation  considered  here  (Equations  (13) )  and  permits 
determination  of  the  range  dependence  and  beam  broadening  of  EM-waves  traveling 
through  vegetation.  Study  of  this  equation,  therefore,  appears  well  justified. 
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although  an  investigation  of  depolarization  effects  would  require  consideration 
of  the  three  remaining  equations  as  well. 

The  conditions  under  which  decoupling  of  the  equation  for  the  first  Stokes 
parameter  occurs  can  be  expressed  in  terms  of  the  co-polarized  and  cross-polar¬ 
ized  field  strength  components  of  scattered  wave  trains  and  their  correlation. 
Though  conceptually  not  difficult,  a  detailed  discussion  of  these  conditions  is 
rather  tedious  and  goes  beyond  the  scope  of  this  report.  In  physical  terms,  the 
conditions  can  be  stated  simply  by  considering  a  wavetrain  incident  from  direction 
T'  and  scattered  into  direction  s.  Decoupling  will  occur  when  the  total  intensity 
of  the  scattered  wave  train  will  depend  only  on  the  total  intensity  of  the  inci¬ 
dent  wave  train  but  (when  averaged  over  many  s'  s  scatter  events  with  both  s' 
and  s'  fixed)  is  independent  of  the  polarization  properties  of  the  incident  radi¬ 
ation.  It  can  be  assumed  that  this  condition  is  satisfied  for  mm-waves  propagat¬ 
ing  in  a  forest  environment. 

3 •  Sol uti on  of  the  Problem. 

We  write  with  Equation  (7) 


P(y)  =  otq(y)  +  (1  -  a) 

(14a) 

where 

2 

.(  X-) 

q(y)  =  (~~)ZQ  Ays  ,  Ay_«ir 

Ay.  s 

(14b) 

is  the  forward  lobe  of  the  phase  function  p  and  (l-<x)  is  the  isotropic  background. 
To  solve  Equations  (13),  we  split  the  diffuse  intensity  Id  into  two  parts 

Id(Z;0 )  =  IjlZ;  04)  +  I2(Z;0)  (15) 


where  I ^  is  determined  primarily  by  the  forward  lobe  q(y)  of  the  scatter  function 

and  I ^  is  the  remainder  of  Id  which  is  primarily  determined  by  the  isotropic 

background.  Thus,  Equation  (13c)  is  separated  into  two  equations: 

I  a  ^  11 

jr  i  r  r 

cos  0-52  +  («A  +  °s)  I J  =  4“ay  J  q  (  y)  Ii  U;0V|>‘)  sino  'do  'd<J> ' 

O  0 


A  +  qS 
Cos°p 


q(Yp) 


(16a) 
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2  it  tt 


cosO  |y2  +  (aA  +  ci$)  I2  =  ff  q(y)  I2  UiO")  sin  G'dO'd^ 


o  o 


2^  7i 

f 

J 
0  0 


(16b) 


+  4?  ^  "  a)  fj  ( 1  j  ( 2 ;© '  )  +  1 2 ( Z ; o  " )}  s i n  o" do  ' dt|> " 


+  ¥(1'a)  S,)e 


qA  *  % 

COS0.. 


for  7  >  O 


Equations (16) can  be  interpreted  on  physical  grounds  and  can  be  simplified  by 
introducing  the  following  approximations: 


(a)  Since  the  specific  intensity  1^  is  generated  by  forward  scattering 
(only  the  forward  lobe  q  of  tne  scatter  characteristic  appears  in 
Equation  (16a)), it  can  be  expected  that  Ij  will  be  characterized  by 
a  strong  and  relatively  narrow  peak  about  the  incidence  direction 

of  the  primary  wave,  while  1^  will  be  very  small  in  other  directions. 
Hence,  the  factor  cos  6  in  the  first  term  on  the  left  side  of 
Equation  (16a)  can  be  replaced  by  cos  9  which  simplifies  solution 
of  this  equation  considerably. 

(b)  The  specific  intensity  l?  is  primarily  determined  by  the  isotropic 
background  of  the  scatter  characteristic  p  and  will  have  a  broad 
directional  spectrum.  Hence,  in  evaluating  the  first  integral  on 
the  right  side  of  Equation  (16b),  the  function  q  (y ) ,  which  is 
characterized  by  a  narrow  peak,  can  be  treated  as  a  6 -function. 

With  Equations  (4)  and  (8b),  we  find: 


2rr  it 


2tt  tt 


-4-  f f  q(Y)  I2  (Z;0\sv)  sin  O'dO'd*-'  *I2(Z;0,<|>) 


(10a) 


0  0 


(c)  In  Appendix  A,  it  is  shown  that  with  approximation  (a) 


ff  h  <ZJ  sinOdOd,},  =  Sp  {exp  (-  “ “^Z)  -  exp(-  ~^-Z) } 


0  0 


(18b) 

which  simplifies  the  second  integral  on  the  right  side  of  Equation  (16b). 
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With  these  approximations  Equations  (16a)  and  (16b)  reduce  to 


2b  +1 


O  J,  .  | 

lJp  Dr 


-  aW 

1  "  4ir  J 
o 


q(y)  I 'A')  dli -<4' 


+  4?  3pe  l‘P 


r  * 

a1?  ,  ,  -  w 
3x  d  d 


p  .  j  -  v 
J  i2(T;i'')  jh'  < sp<  “p 

-1 


where 


cos  Y  =  cos  (<{>-<j)>)[(l-p2)(lV2)l  2  +  up" 


COS  Yp  =  COS 


p)[(l-p2)(l-pp)]  2  +  ppp 


(19a) 


( 19b) 


We  have  introduced  here  the  following  normalizations  which  are  in  accordance  with 
common  usage: 

=  co^  ^ 

L.l  “  V 


(20a) 


pp  =  cos  op 


V  3  (oA  +  o$)Z  -  ot Z  (optical  density) 
C  =  (oA  +  (1-a)  os)Z  =  O^Z  =  x  (1-aW) 

(Albedo) 


(20b) 


W  =  —+  — 
°A  as 


“-^TiTTaTS;  -  ■ 


(20c) 


1-aW 


Hath  1^  and  I-,  are  assumed  to  satisfy  the  boundary  conditions  (13a)  and  (13b), 
i .e. , 

x  ~  0  0<y<l 

I  A(  l  ;  M»d>)  =  0  for  -j  (21a) 
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r  =  0  0  <  |i  <  1 

I2(t;u)  =  0  for  (21b) 

T  -►  cu 

Note  that  Equation  (19b)  is  the  well-known  transport  equation  for  isotropic 
scattering.  The  only  difference  is  that  the  scatter  cross  section  is  reduced  by 
multiplication  with  (1-a)  (only  the  fraction  1-a  of  the  total  power  scattered  in 
any  scatter  event  is  transferred  into  the  isotropic  background).  The  solution  of 
the  isotropic  transport  equation  depends  only  on  the  coordinates  x  and  p  (i.e.,  Z 
and  o)  but  is  independent  of  4>.  Hence,  ^  as  indicated  already  in 

Equation  (15).  In  contrast,  Ij  will,  in  general,  depend  on  all  three  coordinates 
t,  p,  and  <J>.  The  only  exception  is  the  case  of  normal  incidence  of  the  primary 
plane  wave  (0p  =  0)  where  full  symmetry  about  the  Z -direction  is  maintained  and  1^ 
is  independent  of  <J>. 

Solution  of  Equation  for  1^ 

Equation  (19a)  is  solved  by  the  method  of  undetermined  coefficients,  that  is, 
by  writing  Ij  in  the  form: 

(i;p,<f>)  ^  (bjt  +  +  b^r"  + - )e~  p  (22) 


where  the  are  unknown  coefficients.  Note  that  representation  (22)  already  i m* - 
isfies  the  boundary  condition  (21a),  that  1^  should  vanish  at  x  =  0  and  >».  The 
in  general,  are  functions  of  p  and  ({>. 


A  recursion  relation  for  these  coefficients  is  obtained  by  using  Equation 
(22)  in  Equation  (19a)  and  observing  that  the  latter  equation  must  be  satisfied 
for  any  r  >  0: 

Mu,*)  -  i  “r“  'I  Or.)  (23a) 


•>, (>•.»)  ■  ?fr  Hr  <1  <V 

2n  +1 

■  it  f  in  /  /  »<r>  d"'d'1, 

P  0  -1 


(23b) 


m  =  2,3,  -  -  - 
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With  the  assumed  Gaussian  forward  lobe  q(y)>  Equation  (14b),  the  recursion  relation 
can  be  solved  in  good  approximation  in  closed  form.  Details  are  shown  in  Appendix 
B ;  the  result  is 


i  i.  _  'j>  1  /aW.iii 

bH)(u,v)  -  4Tt  (  -) 


%,{y 


111  =  1,2  -  - 


(24a) 


where 


q"|(V  1  nv\v 


.  Lm 

in' Ay. 


(24b) 


The  approximation  is  valid  for  ai$  <<■  n,  the  case  considered  here.  Note  that  the 
functions  q()  have  the  same  general  form  as  the  forward  lobe  q  of  the  phase  function 
but  the  beamwidth  &y$  is  broadened  to  V niAy s . 

By  combining  Equations  (22)  and  (24)  we  obtain  for  1^ 


M 


—  1  ,ctW  ,  m  ,  . 

p  irf1  W 


I1(li)l,<l>)  *  e  r  ml  ,-p 

with  cos  y  =  cos  (<(,-<1,  )[{l-p2)(l-p2)]^  +  vmj 


(25) 


M  is  a  suitably  chosen  upper  limit.  The  general  behavior  of  lj.  as  expressed  by 
Equation  (25),  is  in  accordance  with  expectation.  For  small  distances  into  the 
scatter  medium,  the  first  term  of  the  sum  dominates,  which  means  that  Ij  increases 
linearily  with  t.and  its  directional  spectrum  is  narrow.  (1^  is  generated  by  for 
ward  scattering  of  the  coherent  field  component  in  accordance  with  the  scatter 
function  d(Yp))- 


As  i  increases,  Ij  reaches  a  maximum  and  then  decreases  again  because  of  the 
exponential  factor  exp  (-i/|i  ).  Simultaneously,  higher  order  terms  in  the  sum  on 
the  right  side  of  Equation  (25)  become  significant.,  which  has  the  effect  that  (a) 
the  attenuation  rate  is  slowed  down  and  (b)  beam  broadening  occurs. 

The  reduction  of  the  attenuation  rate  can  be  determined  explicitly  by  con¬ 
sidering  the  forward  direction  where  Yp  =  0.  For  this  direction  the  summation  can 
be  performed  in  closed  form.  We  have 
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s 

Ij  ( i;p,«1>)  r  ^  (  l  e 


4  if 


i  _  I 

"p-  o  *  p  )  qM^vp^ 


-  <■  M 

+  e  Mp  \  p^V11 

in- 1  m!  v|i  '  1  1 


(26) 


P 

As  before,  it  is  assumed  that  M  is  sufficiently  large  so  that  the  value  of  1^  no 
longer  changes  when  M  is  increased.  The  minimum  required  M  will  increase  with 


Solution  of  liquation  for  1., 

Lquation  (19b)  for  1^  is  the  well-known  transport  equation  for  isotropic 
scattering.  The  scatter  cross  section  per  volume  which  appears  in  this  equation 
is  determined  by  the  isotropic  background  of  the  phase  function  and  is  given  by 
(l-n)o  .  This  modification  does  not  change  the  mathematical  character  of  the 
equation. 

A  number  of  solution  methods  for  the  isotropic  transport  equation  are  avail¬ 
able  in  the  literature;  see,  for  example,  References  [3-6].  Numerous  additional 
references  can  be  found  in  [3] .  In  the  present  paper  we  use  a  moment  method  ap¬ 
proach  [6]  which  expands  the  p-dependence  of  into  a  series  of  overlapping  tri¬ 
angular  basis  functions.  These  functions,  Fn(x),  are  defined  as  indicated  in  fig.  6. 
The  range  -1  >  p  +  i  is  subdivided  into  N  intervals 


Vl  -  11  -  "n  n  =  1,2  ....  N 


where  N  is  assumed  to  be  odd  and  pn  is  defined  by 


p  =  cos  (1-  j^)n  =  -  cos  ( j^n )  ,  n  =  0,1,2  ....  N 


(27a) 

(27b) 


I  hus; 


for  -  1  «-  i.  <  n 


n- 1 


V'Vi 

F"<1,)=  '  'V.- 

pn+l"Mn 


ii -1 


n  •-  P. 


|1  <  M  <  11  , 

-  n-1 


(28a) 


nil 


n  =  1,2  ....  N-1 
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The  "boundary"  functions  F  (x)  and  FN(x)  have  the  shape  of  half-triangles  and  are 
defined  by 


F0<n  > 


for  -1  <  p  <  p. 


p 


(28b) 


Fm(m)  = 


-1  <  p  <  p. 


Vl  ,J  -  +1 


(28c) 


Expansion  of  the  specific  intensity  Ij  into  the  functions  Fn(p)  yields 

■  l  an  (f)  Fn(„) 


where  the  expar  ion  coefficients 

an(f)  =  I2  (t;pn)  0,1,2  ....  N  (29a) 

are  the  values  of  I„  at  p  =  p  .  In  other  words,  Equation  (29a)  replaces  I~  by  a 

C  D  *- 

piecewise  linear  approximation  in  p,  which  is  exact  at  the  points  p  =  Pn. 

With  Equation  (29a),  the  transport  equation  (19b)  reduces  to  a  system  of 
coupled  differential  equations  of  the  first  order  in  x 

*  *„  *  I  •  j|„0'Vk>  *  “r  V  (30a) 

n  =  0,1,  ....  N 


where 


zW+i "  *Vi>=  sin  ^  si"  (T} 


J'  Fk(p)dp 


for  k  -  1,2,  .... N-l  (30b) 


\  (1+Pi)  =  \  (  )  =  sin2  (^j) 


S 


We  have  assumed,  furthermore,  that  the  direction-cosine  up  of  the  incident 
plane  wave  coincides  with  one  of  the  u  values  (n  =  j);  for  normal  incidence,  in 
particular,  we  have  j  =  N.  Equations  (30),  a  linear  system  of  coupled  differential 
equations  of  the  first  order,  can  be  solved  by  assuming  an  exponential  T-dependence 
of  a  .  It  is  straightfoward  to  show  that  the  general  solution  takes  the  form: 


an  ™  =  £ 


_  _X.. 

e 

P. 

J 


N 

+  L 
k=o 


Ak 

— * —  e 


n  =  0,1,2 


°jn  =1  0 


for  n  =  j 
n  7  j 


The  first  term  within  the  parentheses  is  a  particular  solution  of  the  inhomo¬ 
geneous  equation,  the  sum  term  is  the  general  solution  of  the  homogeneous  equation. 
The  coefficients  An  are  unknown  amplitudes  which  have  to  be  determined  from  the 
boundary  conditions.  The  attenuation  coefficients  Sk  are  solutions  of  the  char¬ 
acteristic  equation 


-  T 


=  l 


which  has  to  be  solved  numerically.  Useful  insight  into  the  general  behavior  of 
the  spectrum  can  be  gained  by  plotting  the  left  side  of  Equation  (32)  as  a 
function  of  S.  From  Equations  (27b)  and  (30b),  it  can  be  seen  that 


p  =  p 
N-n  r  n 


and  the  left  side  of  Equation  (32),  L(s),  is  symmetric  in  S.  The  same  is  true,  of 
course,  for  the  right  side  R(S)  =  1. 

In  Fig.  6,  L(S)  is  sketched  for  N  =  9.  The  function  has  N  +  1  poles, 
located  at  the  points  S  =  pn-  In  the  regions  between  adjacent  poles  L(S)  varies 
from  -  °°  to  +».  (The  only  except  ion  is  the  interval  between  the  two  poles  closest 
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to  S  =  0.)  Therefore,  L(S)  crosses  R(S)  =  1  N-!  times  between  -1  <  S  <  1.  Further¬ 
more,  for  S  +  +  oo,  We  have* 


,  N 

L(S)  -  W  s  Pn  =  W 
2  n=o  n 


Since  W  <  1,  L(S)  intersects  the  line  R(S)  =  +1  twice  for  |S|  >  1.  Fig.  6 
indicates  the  N+l  solutions  of  Equation  (32)  which  are  distributed  symmetri¬ 
cally  about  the  point  S  =  0.  Counting  these  solutions  in  order  of  ascending  value, 
the  first  N+J[  eigenvalues,  S  . SM  ..  ara  negative  and  would  lead  to  a  specific 

g  O  li  1  t 


intensity  which  increases  exponentially  with  x-  The  boundary  condition  at  £  -*■  00 

eliminates  solutions  of  this  type.  The  remaining  N+l  solutions,  S^+j . S^,  are 

positive  and  lead  to  physically  meaningful  results.2  T~ 


Fig.  6  shows  that  all  Sk  have  magnitudes  smaller  than  unity  except  for  the 
largest  eigenvalue  which  is  >1  and  provides  the  lowest  attenuation  rate.  The 
associated  term  in  Equation  (31)  is  of  particular  interest  since  it  determines  the 
asymptotic  behavior  of  at  large  distances  into  the  scatter  medium. 

The  number  of  eigenvalues  Sk  within  the  range  0  <  S  <  1  increases  with  N, 
the  order  of  the  approximation.  For  N -*■<»,  the  Sk  form  a  continuous  spectrum 
covering  the  entire  range  from  0  to  1.  The  eigenvalue  S^,  however,  remains  dis¬ 
crete  and  for  N  -+■  «°  approaches  a  limiting  value  S,  >  1  which  will  be  further  dis- 

a 

cussed  below.  Numerical  evaluations  have  shown  that  if  N  is  reasonably  large  (21 
or  greater),  Sw  approximates  S  closely  and  does  not  vary  significantly  as  N  is 

Pi  3 

further  increased. 

As  stated  before,  the  boundary  condition  at  t*-°°  eliminates  all  expansion 
terms  in  Equation  (31)  associated  with  negative  eigenvalues  Sk>  Me  have 

A.  =  0  for  k  =  0,1  -•  Nj^l  (33a) 

k  2 


The  remaining  amplitudes  . A^  are  determined  by  the  boundary  condition  at 

~T 

•t  =  0.  By  applying  Equation  (21b)  to  Equations  (29)  and  (31)  the  following  system 


N 

*From  Fig.  5,  it  is  easy  to  see  that  ^  Fn(p)  =  *  f°r  ^1^5  +!• 


n  r  n 

Hence  2  =  /  2  F ( jj ) djj  -  2. 

n=o  J  n=o 
-1 
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of  linear  equations  is  obtained  for  these  amplitudes: 

N 


Un  -a-  ■  ito 

K  2  1  lJn  P. 


,  N+l  .. 

for  n  =  -rr~  ""  N 


(33b) 


This  system  has  to  be  solved  by  computer.  With  the  A^  and  determined,  the 
solution  for  takes  the  final  form 


I  o  ( T «  p )  =  J>  { 
2ir 


Pi 


-  e 


Fjtu) 

"PJ 


+  Jn±1  [  flke 

K”  o 


•  ,5C- 

k  F(m) 


1-  un 
S". 


-  1  > 


(34) 


for  x  >  0,  -1  <  y  <  +1 


Pj  -  cos  ep  =  yp 


It  is  instructive  to  examine  the  behavior  of  I2  at  small  and  very  large  x- 
For  small  T ,  I2  can  be  approximated  by 


I2  (t;p) 


for  0  <  p  =  cos  6  <  1 


-1  <_  p  =  cos  e  <  0 


(35) 


This  so-called  first  order  multiple  scattering  approximation  is  obtained  if 

Equation  (19b)  is  solved  with  the  integral  term  on  the  right  side  suppressed  [3]. 

Equation  (35)  shows  that  is  zero  at  the  air-forest  interface  for  0  <  6  <  tt,  as 

c  ~  2 
required  by  boundary  condition  (21b).  As  x  is  increased,  I2  at  first  increases 

linearly  with  distance  (I2  is  generated  by  scattering  of  the  coherent  field  com¬ 
ponent  Iri  in  this  region)  but  then  reaches  a  maximum  and  decreases  again  because 
of  absorption  in  the  scatter  medium.  At  very  large  distances,  I2  is  determined 
by  the  asymptotic  solution  [3]: 
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The  quantities  S  .  and  f,  are  shown  in  Figs.  7  and  8  as  functions  cf  W.  The 

d  d  A 

asymptotic  amplitude  A^  depends  on  the  boundary  conditions  and  is  plotted  versus  W 
in  Fig.  9  for  the  conditions  shown  in  Equation  (21b).  The  figures  show  that  in  a 
strongly  scattering  medium  with  ft  close  to  unity,  the  asymptotic  attenuation  rate 
1/S  is  low,  the  amplitude  A  is  substantial,  and  the  directional  spectrum  f  is 

da  a 

wide.  Experimental  results  indicate  that  the  forest  environment  has  these  char¬ 
acteristics  and  is  a  strongly  scattering  medium  at  mm-wave  lengths.  In  Appendix  C 
it  is  shown  that  the  expansion  term  in  Equation  (34)  associated  with  the  largest 
eigenvalue  converges  towards  Equations  (36)  as  N  •+•  <*>. 

A  final  remark  is  concerned  with  backscattering  out  of  the  forest  medium 
into  the  air  halfspace.  While  J  .  is  zero  and  1^  is  negligible  in  the  range  of 
the  backscatter  directions  -1  <  u  =  cos  e  <  0,  I.,  is  different  from  zero  in  this 
range.  3ackscatter  into  the  air  halfspace  is  determined  by  the  u-spectrum  of  I., 
at  x  =  0.  With  Equation  (35),  we  have 

r 

l„(0;u)  =  _p_  ft  for  -1  <  p  =  cos  0  <  0 

4ir  1-p 

Apparently,  side  scatter  (o  =  90°)  is  moderately  stronger  than  direct  back¬ 
scatter  (o  =  180°).  Note  that  the  backscatter  radiation  is  generated  by  the 
(isotropic)  scattering  within  the  forest  medium;  but  it  is  not  caused  by  reflection 
of  the  incident  energy  at  the  air-to-forest  interface  (which  is  insignificant).  It 
is  a  volume  rather  than  a  surface  effect. 
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4 .  Power  Received  by  Highly  Directive  Antenna 

The  present  theory  assumes  that  the  forest  halfspace  is  illuminated  by  a 
plane  wave  incident  from  the  air  halfspace.  This  situation  can  be  approximated  by 
locating  the  transmit  antenna  at  a  sufficiently  large  distance  from  the  edge  of 
the  forest.  We  assume,  furthermore,  that  the  receive  antenna,  which  is  located  in¬ 
side  the  vegetated  area,  is  a  highly  directive  antenna  of  narrow  beamwidth.  This 
is  in  accordance  with  the  experiments  to  which  the  theory  will  be  compared  in  the 
next  section.  We  characterize  the  receiving  antenna  by  its  power  radiation  pattern 
Gr  ( Yr )  where  yR  is  the  angle  which  the  direction  of  observation  6,*  includes  with 
the  pointing  direction  of  the  antenna  axis  0R,*R  (main  beam  direction);  see  Fig. 

10.  Evidently 

cos  YR  =  cos  ( ^ q )  sin  0  sin  0R  +  cos  0  cos  eR  (37a) 

We  normalize  GR  such  that 

2ir  it 

'  ^ GR(yR)  sin  0  d  0  d  J>  =  4n  (37b) 

o  o 

The  directivity  gain  of  the  antenna  is  then  given  by  GR(0).  The  power  PR 

received  bv  the  antenna  is  determined  from  the  well-known  relation  A  -  *o  G  for 

ett  W 

the  effective  receiving  cross  section  of  a  perfectly  matched  antenna  of  directivity 
G.  Thus 

2n  +1 
\  £  r  r 

PR(X;  V  V  =  j  J  Gr(Yr)  I  (T;  p,*)  dud*  (33) 

where  X  is  the  wavelength.  We  have  used  here  the  normalized  coordinates 
t  =  (o^  +  0^)  Z,  pR  =  cos  0R  and  p  =  cos  o.  PR,  of  course,  depends  on  the  location 
and  orientation  of  the  antenna. 

From  the  analytical  viewpoint.,  it  is  convenient  to  assume  a  Gaussian 
radiation  pattern 


w "  (4}  e 


,<Y1  )2 
W 


(39) 


where  AyD  is  the  beamwidth  of  the  antenna.  The  receive  antennas  used 
K 


in  the  experiments  had  a  very  small  beamwidth  of  1.2  and  we  shall  assume  in  the 
numerical  evaluation  that  the  width  of  the  forward  scatter  lobe  of  the  phase 
function  of  the  medium,  Ays,  exceeds  AyR  by  a  factor  of  3-10.  (The  aperture  size 
of  the  test  antennas  was  significantly  larger  than  the  cross  section  of  most  scat¬ 
ter  objects  in  a  forest  environment  with  the  exception,  perhaps  of  tree  trunks.) 

Using  Equations  (12),  (26)  and  (34)  for  1^  and  I^,  respectively,  we  obtain 
for  the  received  power 


/^RPv2  T 

exp  {-  (t--)  -  -f-  } 

AyR  yp 


+  “  {  [e  "  e  1  <!„<%>> 

j_  M  i  u  m  (40a) 

+  e  Up  m-l  tqm(YPR)‘qr4YPR)l1 

2  ~  7~  N  $.  Pn^R^ 

U|>  +  'tV  V  1» 

Pj  ^  ^o'-Sk 


where  g.  (ypR)  -  -7-— 7 
\r  *  "“yS 


%R  +  n,AyS 


2  2  ' 

COS  Ypp  =  cos  (i}>p-ti>R)  l  ( l'Pp)  ( 1~PR)  1  ^  +  PpPR 


(40b) 


(40c) 


We  have  normalized  the  received  power  by  division  with  the  maximum  power 
P,1!lu  which  is  received  if  the  antenna  is  placed  at  the  forest  boundary  t  =  0  and 

IllaX 

is  aligned  with  the  incident  radiation,  i.e.,  0R  =  Op,  -  4>p.  Evidently 

2  2 

Pmax  =  V0;  "P’V  =  47  tiR(0)Sp  =  47  (jjH  $p 

yK 

The  first  term  on  the  right  side  of  Equation  (40a)  is  the  contribution  of  1(,. 
to  the  received  power,  the  second  and  third  terns  are  those  due  to  Ij  and  I^.  The 
integrals  occuring  in  the  contribution  of  I.  have  been  evaluated  by  using  Equation 


(B.6),  and  the  integrals  in  the  term  due  to  1^  were  solved  by  regarding  GR  as  a 
6-function.  This  is  justified  since  is  independent  of  <h  and  changes  very  slowly 
withn  as  compared  to  uR.  The  normalized  variables  are  defined  by  Equations  (20). 

At  small  distances  into  the  random  medium,  for  t  in  the  order  of  3  or  smaller, 
the  coherent  field  component  dominates  and  l’R/i’|)wx  is  in  good  approximation  given 
by  the  first  term  of  Equation  (40): 


P 

max 


-  (>)2 


M  Y  |>r 


) 


(41) 


In  this  region,  the  angular  dependence  of  PR  is  detev'mined  by  the  radiation 
pattern  of  the  receiving  antennas  and  the  attenuation  rate  t/Z  -  +  o<.  is  com¬ 
paratively  high.  Both  attenuation  and  scatter  reduce  the  power  of  coherent  com¬ 
ponent. 


In  the  range  of  intermediate  distances  (r  in  the  order  of  10),  the  forward 
scattered  component  1^  dominates  and  PR  is  primarily  determined  by  the  second  term 
in  Equation  (40).  In  this  region,  the  beamwidth  gradually  broaoens  as  higher  and 
higher  terms  in  the  sum  over  M  become  significant  while  the  attenuation  rate  de- 
c reuses  from  i/Z  =  oA  t-  os  to  c/7.  =  oA  +  (l-«)  o s>  I,  vs  attenuated  due  to  ab¬ 
sorption  and  scattering,  but  scattering  into  the  isotropic  background  only. 


At  large  i ,  the  received  power  is  primarily  determined  by  the  contribution 
of  1 2 *  i.e.,  by  the  last  term  in  Equation  (40).  Asymptotically,  PR  is  with  Equa 
Lion  (36a)  given  by 


max 


2Aa 

Gj^O) 


(Ur) 


fo> 


(42) 


Afl,  fa(„)  and  i/sa  are  shown  in  Figs.  /  to  9  as  functions  of  W.  For  large 
W  1,  the  practically  interesting  case,  the  directional  spectrum  is  broad  and 
1/S  is  small  so  that  the  asymptotic  attenuation  rate  t/(Z-$J  -  (on  (  (l--u)oc)/3, 

t*  a  n  o  a 

is  significantly  reduced  from  the  initial  value  of  oA  +  Oy  The  intensity  1^  de¬ 
creases  essentially  due  to  absorption  only  (though  over  on  extended  path  length 
because  of  multiple  scattering);  but  it  scatters  into  itself. 


In  summary:  While  a  mm  wave  beam  transmitted  through  a  forest  will  he  highly 
directive  and  strongly  attenuated  at  short  distances,  't  will  transform  into  a 


radiation  field  of  large  beaniwidth  attenuated  at  a  much  reduced  rate  at  large  dis¬ 
tances.  These  theoretical  predictions  are  in  agreement,  with  expectation  and  have 
been  confirmed  by  the  experiments. 

A  similar  range  dependence  in  woods  and  forests,  i.e.,  a  transition  from  a 
high  attenuation  rate  at  short  distances  to  a  much  lower  rate  at  large  distances 
has  also  been  observed  at  much  lower  frequencies,  i.e.,  for  propagation  of  VMF 
signals  through  dense  vegetation  (jungles).  The  underlying  mechanism,  however, 
is  completely  different.  The  change  of  attenuation  rate  of  VHF  signals  is  attrib¬ 
uted  to  the  so-called  up-over-down  propagation  mechanism  where  a  lateral  wave  is 
excited  which  travels  along  the  horizontal  boundary  surface  between  the  forest 
canopy  and  the  air  region  above  it  (Tamir  theory).  Since  this  wave  propagates 
mostly  in  the  air  region,  it  is  much  less  attenuated  than  the  direct-wave  travel- 
1  ing  through  the  forest  medium. 

However,  the  lateral  wave  has  a  substantial  launching  loss.  Hence,  at  short 
distances,  the  power  transfer  between  two  antennas  located  in  the  jungle  is  deter¬ 
mined  by  the  direct  wave  and  the  attenuation  increases  rapidly  with  antenna 
spacing.  U  large  distances,  on  the  other  hand,  the  lateral  wave  will  dominate 
and  the  attenuation  rate  is  significantly  reduced.  Note  that  the  effect  is  deter¬ 
mined  by  the  coherent  field  component  alone  while  the  incoherent  component  is 
negligible  at  wavelengths  in  the  VHF  region.  In  contrast,  in  the  imn-wave  region, 
the  incoherent  component  is  significant  in  forests,  and,  as  explained  above,  will 
dominate  at  large  distances.  The  decreasing  attenuation  rate  is  caused  by  the 
interplay  between  this  component  and  the  more  strongly  attenuated  coherent  compo¬ 
nent.  A  lateral  wave  supported  by  the  canopy-to-air  interface  should  not  be  ex¬ 
pected  at  mm-wave  length.  The  interface  is  extremely  rough  in  the  mm  scale. 
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III.  NUMERICAL  EVALUATION 


A  compute,'  program  for  the  evaluation  of  Equation  (40)  has  been  written. 

While  the  first  and  second  terms  of  this  equation  (which  represent  the  contribution 
of  1^  and  1^)  appear  in  explicit  form,  the  third  term  contains  the  eigenvalues 
and  amplitudes  A^  of  l?,  which  have  to  be  determined  by  numerical  solution  of  the 
characteristic  Equation  (32)  and  the  linear  system  (33b),  respectively.  A  null 
finding  procedure  and  a  matrix  inversion  subroutine  for  the  solution  of  these 
equations  are  included  in  the  program. 

Numerical  results  are  shown  in  Figs.  11  through  20.  Table  1  summarizes  the 
parameter  values  applying  to  the  various  curves.  Throughout  it  is  assumed  that  the 
receiving  antenna  has  the  same  3  dB  beamwidth  of  1.2°  as  the  antenr.a  used  in  the 
experimental  study,  which  means  that* 

AyR  =  0.7°  =  0.012  rad 

Furthermore,  we  assume  that  the  3  dB  width  of  the  forward  scatter  lobe  of  the 
phase  function  of  the  scatter  medium  is  6°,  or 

Ays  =  3.5°  =  0.060  rad 

The  quantity  Q  appearing  in  the  figures  is  the  normalized  received  power  expressed 

in  dB,  i ,e. ,  p 

Q  =  10  lg(7r-  ) 
max 

The  distance  into  the  forest  halfspace  is  expressed  in  terms  of  the  normalized 
vegetation  depth  t  =  (o^  +  a<-)Z  as  measured  from  the  air-to-forest  interface;  t  is 
sometimes  referred  to  as  "optical  distance." 

Discussion 

Figs.  11-13  show  the  range  dependence  of  the  received  power,  i.e.,  Q  =  Q ( t ) 
for  various  values  of  W  and  a.  It  is  assumed  here  that  the  primary  plane  wave  is 
incident  normal  to  the  air-to-forest  interface  (plane  Z  =  0)  and  that  the  receive 
antenna  is  pointed  into  the  direction  of  this  wave;  thus,  6p  =  0R  =  0.  The  curves 
demonstrate  the  change  in  attenuation  rate  from  a  high  value  at  small  r  to  a  lower 
value  at  l„rge  c.  The  effect  is  pronounced,  in  particular,  for  strongly  scatter¬ 
ing  media  with  o-  »a^  where  the  albedo  W  is  close  to  unity.  Increasing  a,  the 

*AyR  is~the  anglir counted^  fromTiTe  antenna  axis  at  which  the  mainlobe  of  the  power 
radiation  pattern  of  the  antenna  is  down  by  a  factor  1/e;  see  Equation  (39). 

This  angle  is  related  to  the  3  dB  beamwidth  of  the  antenna  by  A  =0.6  A  .. 

A  <.  is  similarily  defined. 
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(text  continues  on  page  50) 
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FIGURE  13:  Range  Dependence  of  Received  Power  for  Normal 
Incidence  and  a  =  0.8 
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FIGURE  14b :  Directional  Dependence  of  Received  power  for  Normal 
Incidence  and  a  =  0.5,  W  =  0.75 
(b)  Expanded  s- Scale 
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FIGURE  16a:  Directional  Dependence  of  Received  Power  for  Normal 
Incidence  and  a  =  0.6,  W  =  0.95 
(a)  360°  Scan 
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FIGURE  16b:  Directional  Dependence  of  Received  Power  for  Normal 
Incidence  and  a  =  0.5,  W  =  0.95 
(b)  Expanded  ©-scale 
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Q  (dB) 


Directional  Dependence  of  Received  Power  for  Normal 
Incidence  and  u  s  0.2,  W  =  0.9 


45 


Dependence  of  Received 
e  and  q  -  0.8,  W  =  0.9 
Expanded  e- scale 


ratio  of  the  forward  scattered  power  to  the  total  scattered  power,  has  a  similar 
effect,  which  is  easy  to  understand  since  the  attenuation  rates  of  the  field  inten¬ 
sity  components  1^  and  l£,  which  dominate  at  large  r,  are  proportional  to  (1-aW) 
in  this  region. 

Typically,  for  W  =  0.75  and  a  =  0.8,  the  coherent  field  component  I  .  oro- 
vides  the  dominant  contribution  to  the  received  power  in  the  range  0  <  t  <  7, 
while  is  the  dominant  component  in  the  range  7  <  r  <  50,  and  in  the  range 
i  >  50.  If  W  is  increased,  i.e.,  for  W  =  0.95,  a  =  0.8,  the  changeover  occurs 
earlier,  i.e.,  at  t  =  5  rather  than  7,  and  30  rather  than  50. 

The  curves  Q  =  Q(0R)  of  Figs.  14  to  18  show  the  dependence  of  the  received 
power  on  the  pointing  direction  of  the  receiving  antenna.  It  is  assumed  here  that 
the  antenna  is  placed  at  various  vegetation  depths  r  where  it  performs  a  full  360° 
scan.  As  before,  it  is  assumed  that  the  primary  plane  wave  is  normally  incident 
upon  the  air-to-forest  interface.  In  this  case,  Q  depends  on  0R  only  but  is  inde¬ 
pendent  of  cj>R.  Figs.  14b,  16b  and  18b  show  Q  (eR)  in  an  expanded  eR-scale. 

For  t  =  0,  when  the  receiving  antenna  is  positioned  in  the  interface  plane, 
the  curves  Q  (9R)  show  the  Gaussian  radiation  pattern  of  the  receiving  antenna. 

(The  incident  field,  in  this  case,  is  an  unperturbed  plane  wave.)  This  holds  for 
the  two  forward  quadrants  0  <  0R  <_  90°  and  270°  <  0R  <  360°.  In  the  range  of  back 
ward  directions,  90°  1  0  £  270°,  the  antenna  receives  the  backscatter  radiation 
emerging  from  the  forest.  The  figures  show  that  this  backscatter  radiation  has  a 
rather  uniform  directional  distribution  with  sidescatter  moderately  higher  than 
direct  backscatter,  i.e., 

Q  (£)  =  Q(-y-)  ~  Q(tt)  for  T  =  0 

As  T  is  increased,  beam  broadening  occurs.  For  small  r  =  0.1  and  0.3,  the 

coherent  component  remains  dominant,  as  indicated  by  the  narrow  spike  about  0R  =  0. 

But  the  incoherent  component  produces  a  broad  tail  consisting  of  a  peak  of  moderate 

height  and  width  extending  to  about  0R  =  10°,  and  a  very  broad  back  scatter  spectrum. 

The  peak  is  determined  by  I ^  and  the  background  by  I^.  Ij  depends  on  A^,.,  the 

width  of  the  forward  lobe  of  the  phase  function;  if  A  ^  would  be  increased  from  3.5°, 

0 

the  value  assumed  here,  to  7  ,  the  peak  corresponding  to  Ij  would  decrease  in 

height  by  6  dB  and  increase  in  width  to  about  0  =  20°.  For  'r  =  3,  the  coherent  com- 

R 

ponent  is  still  clearly  visible,  but  for  t  =  10  it  has  all  but  disappeared  within 
the  incoherent  component.  At  t  =  30,  (Fig.  18),  the  incoherent  component  has 
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taken  over  completely,  the  signal  level  has  dropped  significantly,  and  the  beam  has 
broadened  out  substantially. 

For  2°  <  0R  <  90°,  i.e.,  in  the  range  of  forward  directions  away  from  the 
spike  of  the  coherent  component,  the  received  power  level  Q  first  increases  with  t, 
up  to  about  t  =  3,  and  then  decreases  again.  This  is  explained  by  the  fact  that 
the  incoherent  component,  which  determines  the  received  power  in  this  6-region,  is 
generated  by  the  scattering  of  the  coherent  component  in  the  range  of  small  x,  but, 
after  reaching  a  maximum,  decreases  again  because  of  absorption.  The  maximum 
occurs  near  t  =  x(l-aW)  =  1  which  roughly  corresponds  to  x  =  3  for  most  of  the  a 
and  W  values  considered  here. 

The  receiving  antenna  assumed  here  is  an  idealized  antenna  radiating  a  narrow 
mainbeam  but  no  sidelobes.  Practical  antennas  produce  sidelobes  and  the  measured 
directional  spectrum  Q  (eR)  will  be  influenced  by  the  sidelobe  structure  of  the 
test  antenna.  It  is  likely  that  such  distortions  will  obscure  some  of  the  details 
of  the  directional  spectrum  of  1^  and  in  particular,  those  at  a  level  of  more 
than  30  or  40  dB  below  the  main  peak. 

Figs.  19  and  20  show  the  range  dependence  and  directional  spectrum  of  the  re¬ 
ceived  power  for  the  case  that  the  primary  wave  is  incident  under  an  oblique  angle, 
i.e.,  6p  =  45°,  4>p  =  0.  In  calculating  the  range  dependence  curves  of  Fig.  19,  it 
was  assumed  that  the  receiving  antenna  is  pointed  into  the  direction  of  the  inci¬ 
dent  wave  such  that  0R  =  45° ,  4>R  =  0.*  The  figures  show  the  same  general  behavior 
of  Q,  the  normalized  received  power,  as  those  for  normal  incidence,  with  *he  ex¬ 
ception  that  in  Fig.  20  the  maximum  of  the  received  power  occurs  at  0R  "  45°  rather 
than  0°.  A  second  difference  may  be  seen  by  comparison  of  Figs.  19  and  12,  which 
show  the  range  dependence  of  Q,  at  oblique  and  normal  incidence,  for  the  same 
values  of  a  and  W.  For  W  =  0.5,  when  absorption  is  strong,  the  attenuation  rate 
for  oblique  incidence  is  significantly  above  that  for  normal  incidence.  This  is 
due  to  the  fact  that  in  the  case  of  ohlique  incidence  the  wave  trains  received  at 
a  given  point  x  in  the  medium  have  traveled  a  path  which,  in  effect,  is  by  a  factor 
(cos  Op)'1  longer  than  in  the  case  of  normal  incidence.  Hence,  the  attenuation  is 
that  accumulated  over  a  longer  path.  (Note  that  the  receive  antenna  is  a  narrow 
beam  antenna  pointed  into  the  direction  6R  =  0R  so  that  it  receives  primarily  for¬ 
ward  scattered  wave  trains.)  The  effect  is  clearly  visible  in  the  x-range  where 
the  "directive"  field  components  and  1^  dominate,  but  it  levels  off  in  the 

*The  curves  in  Fig.  19"  shoul d  start  at  Q  =  0  dB  for  x  =  0.  The  deviation  is  due  to 
a  computational  quantization  error. 


range  of  very  large  x  where  1 2  is  the  strongest  field  component.  1 2  propagates  by 
omnidirectional  scattering  and  its  magnitude  at  any  point  t  is  primarily  determined 
by  the  constituent  wave  trains  having  traveled  along  the  shortest  paths,  i.e., 
along  those  directed  normally  into  the  forest  medium.  (But  the  directional  spec¬ 
trum  of  Ij  is  very  broad  so  that  the  received  power  does  not  depend  significantly 
on  antenna  pointing.) 

IV.  EXPERIMENTS 

Systematic  test  series  on  mm-wave  propagation  in  woods  and  forests  were  con¬ 
ducted  by  Violette,  Fspeland,  et.  al.*  Detailed  reports  on  these  studies  have  been 
published  recently  [ 1 ,  2].  The  following  discussion  draws  strongly  from  these 
reports. 

1 •  Description  of  Experiments 

The  experiments  v/ere  conducted  at  three  frequencies  simultaneously,  i.e.,  at 
9.6,  28.8  and  57.6  GHz.  The  first  of  these  frequencies  is  at  X-band;  the  two  remain¬ 
ing  frequencies  are  close  to  the  35  and  60  GHz  bands  which  are  primarily  of  inter¬ 
est  in  tactical  communication.  The  transmit  antennas  have  a  comparatively  broad 
beainwidth  of  10°  and  provide  wide  angle  illumination  of  the  test  area.  The  re¬ 
ceiving  antennas  have  narrow  beamwidth,  i.e.,  4°  at  X-band  and  1.2°  at  the  two  mm- 
wave  frequencies.  In  the  numerical  evaluation  of  the  theory,  the  1.2°  beamwidth 
was  assumed  for  the  receive  antenna. 

The  antennas  are  linearily  polarized  and  in  the  experiments  of  interest  here 
were  operated  in  co-polarization.  The  experiments  have  shown  that  at  mm-wave 
frequencies,  the  propagation  conditions  in  vegetation  are  practically  the  same  for 
vertical  and  horizontal  polarization  so  that  the  direction  of  polarization  is  of  no 
consequence.  This  was  also  assumed  in  the  theory.  There  is,  however,  the  follow¬ 
ing  difference.  The  theory  is  based  on  the  scalar  transport  equation  which  means 
that  the  received  power  is  obtained  as  the  total  power  intercepted  from  both  the 
co-polarized  and  cross-polarized  field  components.  In  the  experiments,  which  util¬ 
ize  linearily  polarized  antennas,  only  the  copolarized  component  contributes  to 
the  received  power.  Cross -polarization  experiments  have  shown,  however,  that  over 
the  distances  considered  here,  mm-wave  beams  do  not  become  fully  depolarized  in 
vegetation  and  that  the  cross-polarized  component  remains  several  dB  below  the  co¬ 
polarized  component.  Hence,  errors  due  to  this  discrepancy  of  theory  and  experi¬ 
ments  should  not  exceed  1-2  dB. _ 

*The  investigators  are  with  NTIA/1TS,  US  Dept,  of  Commerce,  Boulder,  CO.  The  ex¬ 
periments  were  performed  under  contract  with  the  US  Army  Communications- Electron¬ 
ics  Command,  Fort  Monmouth,  NJ. 
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Normally,  the  vegetation  density  in  a  forest  will  change  strongly  with  lo¬ 
cation  and  is  difficult  to  characterize.  To  eliminate  this  uncertainty,  the  experi¬ 
ments  were  conducted  in  a  regularly  planted,  well-groomed  orchard  of  pecan  trees 
(near  Wichita  Falls,  TX)  which  should  approximate  a  statistically  homogeneous  vege¬ 
tation  medium  -  as  far  as  this  is  possible.  The  trees  in  this  orchard  were  of 
roughly  equal  growth  they  were  planted  in  a  regular  square  pattern  spaced  about 
13  m  apart,  and  individual  trees  had  a  height  of  approximately  10  m.  Maximum  vege¬ 
tation  density  occured  in  the  canopy  region  between  4  and  6  m  above  ground.  The 
experiments  were  conducted  both  under  summer  and  winter  conditions,  i.e.,  with 
trees  in  leaf  and  without  leaves.  Fig.  21  shows  the  test  site  in  early  April  and 
late  August  when  the  experiments  were  conducted. 

The  experiments  which  are  of  interest  in  the  context  of  the  present  paper 
can  be  described  with  the  help  of  Fig.  22.  The  black  circles  in  this  figure  indi¬ 
cate  some  of  the  trees.  The  transmit  antenna  was  located  at  or  near  the  point  T ^ 
outside  the  orchard  at  a  distance  of  -300  m  from  the  first  row  of  trees.  The  re¬ 
ceive  antenna  was  located  within  the  orchard  and  was  moved  to  various  locations  R3, 
R^,  R^,  etc.,  as  indicated  in  the  figure.  In  each  case,  the  receive  antenna  was 
positioned  at  the  same  height  as  the  transmit  antenna,  and  the  1 ine-of-sight  follow¬ 
ed  a  path  of  maximum  vegetation  density  through  one  column  of  trees.  The  measure¬ 
ments  were  taken  at  three  different  heights  above  ground,  i.e.,  at  1  m,  4  m,  and 
6  in.  The  4  and  6  m  experiments  are  primarily  of  interest  here.  At  each  receiver 
position  and  height,  the  receive  antenna  was  first  pointed  towards  the  transmit 
antenna  so  that  its  main  beam  direction  coincided  with  the  1 ine-of-sight.  The 
receive  antenna  then  performed  a  scan  in  azimuth  by  +15°  and  in  elevation  by  +10°. 
The  quantity  measured  was  the  received  power  at  constant  transmit  power. 

2 .  Experimental  Results 

Results  applicable  to  this  paper  are  shown  in  Figs.  23  to  26.  In  the  first 

two  figures,  the  range  dependence  of  the  received  power  is  plotted.  The  quantity 
shown  is  the  maximum  of  the  received  power  observed  during  the  azimuth  and  ele¬ 
vation  scans.  The  power  is  normalized  to  0  dB  at  the  edge  of  the  orchard.  Fur- 

_  o 

thermore,  free-space  path  losses  (R  dependence)  and  atmospheric  attenuation  (at 
57.6  GHz)  have  been  subtracted  so  that  the  curves  show  the  vegetation  loss  only. 
Because  of  this  normalization,  comparison  with  the  theoretical  results  based  on 
the  assumption  of  an  incident  plane  wave  seems  appropriate.  The  abscissa  scale 
shows  the  number  of  trees  on  the  1 ine-of-sight;  this  number  multiplied  by  the 
tree  spacing  of  13  in  yields  the  vegetation  depth  in  meters.  The  solid  curves 
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FIGURE  22:  Drawing  of  Test  Area 

The  black  circles  show  some  of  the  trees.  In  he 
experiments  of  interest  here,  the  transmitter  terminal 
was  located  outside  the  orchard  at  point  T3A.  The  re¬ 
ceiver  terminal  was  positioned  inside  the  orchard;  re¬ 
ceiver  locations  are  indicated  by  R^,,  R^,  etc. 

(From  Violette  and  Espeland  [2]. 
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FIGURE  23:  Received  Power  in  dB  {Vegetation  Loss)  as  Function  of  Number 

of  Trees  in  Path. 

The  1 ine-of-sight  (LOS)  followed  a  maximum  vegetation  density  path 
at  a  height  of  4  m  above  ground.  The  trees  are  spaced  *13  m  apart. 

Solid  Lines:  Vegetation  loss  for  trees  in  leaf. 

Dashed  Curves:  Vegetation  loss  for  trees  without  leaves. 

(From  Violette  and  Espeland  [2].) 
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FIGURE  24:  Received  Power  in  dB  (Vegetation  Loss)  as  Function 

of  Number  of  Trees  in  Path 

Height  above  ground:  6  in.  Other  conditions  are  those  of  Figure  23 
(From  Violette  and  Espeland  [2].) 


show  the  simmer  data,  the  dashed  curves  the  winter  data.  Apparently,  leaves  have  a 
strong  effect  on  nm-wave  attenuation.  As  one  should  expect,  the  attenuation  in¬ 
creases  with  frequency  though  the  effect  is  not  dramatic  in  the  nm-wave  region.  It 
is  evident  that  the  experimental  results  confirm  the  theoretical  prediction  that 
the  mm-wave  attenuation  rate  in  woods  decreases  from  a  high  value  at  small  vege¬ 
tation  depths  to  a  much  smaller  value  at  large  depth;  the  transition  between  these 
two  regions  appears  to  be  rather  sharp. 

Figs.  25  and  26  show  the  received  power  as  a  function  of  the  azimuth  and  ele¬ 
vation  scan  angles.  The  light  traces  indicate  the  power  received  in  free  space, 
i.e.,  the  radiation  pattern  of  the  receive  antenna.  The  figures  apply  to  trees 
without  leaves,  but  very  similar  curves  were  recorded  for  trees  with  leaves. 

Again,  good  qual itative  agreement  with  theoretical  predictions  is  observed.  For 
short  distances  into  the  orchard,  i.e.,  at  vegetation  depths  of  one  and  three 
trees,  the  presence  of  the  cohere„t  component  is  evident  by  the  narrow  peak  of  the 
scan  pattern.  At  larger  vegetation  depths,  corresponding  to  propagation  paths 
through  eight  and  eleven  trees,  the  coherent  component  has  disappeared  within  the 
incoherent  component  which  dominates  the  received  power,  resulting  in  a  reduced 
amplitude  and  substantial  beam  broadening. 

The  scintillations  in  the  curves  of  Figs.  25  and  26  are  interference  effects 
(fadings);  the  exact  locations  of  the  maxima  and  minima  should  vary  rapidly  with 
antenna  position  so  that,  in  an  averaged  curve,  these  fluctuations  will  disappear. 
Since  transport  theory  neglects  interference  effects,  the  corresponding  theoretical 
curves  of  Figs.  14  to  18  are  smooth. 

3 .  Comparison  of  Theory  and  Experiments 

Theoretical  and  experimental  results  show  the  same  general  trends  and  are  in 
good  qual  itative  agreement.  One  may  now  try  to  establish  quantitative  agreement 
by  choosing  the  numerical  values  of  the  four  theoretical  parameters  AyS’  V  °S 
and  a,  such  that  a  "best  fit"  of  the  calculated  curves  to  the  measured  curves  is 
achieved.  At  the  same  time,  this  procedure  would  yield  order  of  magnitude  esti¬ 
mates  for  the  four  parameters.  The  procedure,  however,  is  not  without  difficulties 
and,  as  will  become  apparent  from  the  discussion  below,  the  presently  available 
experimental  data  base  is  not  sufficiently  large  to  arrive  at  definite  conclusions. 
It  should  be  understood,  also,  that  the  comparison  of  theory  and  experiments  in¬ 
volves  a  certain  amount  of  conjecturing  since  the  measured  curves,  as  opposed  to 
the  smooth  theoretical  curves,  show  a  significant  amount  of  fluctuation.  Moreover, 
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Directional  Dependence  of  Received  Power  at.  Height 
of  4  m  with  1,  3,  8  and  11  Trees  in  Path. 


Ihe  scan  range  in  azimuth  is  +15°,  in  elevation  ±10°;  the  curves 
apply  to  trees  without  leaves.  The  light  traces  show  the  radiation 
patterns  of  the  receive  antennas. 

(From  Violette  and  fspeland  [2].) 


Elevation  (  °) 

FIGURE  26:  Directional  Dependence  of  Received  Power  at  Height 
of  6  in  with  1,  3,  8  and  11  Trees  in  Path. 

Conditions  are  same  as  in  Figure  25. 


(From  Violette  and  Espeland  [2  |.) 


the  comparison  has  to  be  regarded  with  caution  since  the  theoretical  model  contains 
a  priori  assumptions  and  strong  simplifications  and  is  not  identical  with  the  ex¬ 
perimental  set-up.  In  particular*  the  theoretical  model  applies  to  a  forest  half- 
space  illuminated  by  a  plane  wave,  while  in  actuality  the  incident  beam  has  a 
finite  cross  section  and  the  trees  have  a  finite  height  so  that  losses  occur  due 
to  radiation  into  the  ground  and  into  the  air  region  above  the  forest  canopy;  such 
losses  are  not  accounted  for  in  the  theory. 


In  the  following,  a  critical  assessment  is  attempted  in  how  far  each  of'  the 
four  theoretical  parameters  can  be  determined  from  the  experimental  data.  For  the 
purpose  of  this  assessment,  the  parameters  are  rewritten  in  the  form  Oj  = 
oA  +  W  -  os/(oA  +  os),  and  a. 


(1)  The  width  of  the  forward  lobe  of  the  phase  function  determines  the 
(changing)  beamwidth  and  amplitude  of  the  Ij-component  of  the  incoherent  intensity, 
and  would  have  to  be  determined  from  the  experimental  curves  of  Figs.  25  and  26. 

But  these  curves  are  not  in  a  form  to  facilitate  comparison  with  theoretical  data 
and  no  attempt  is  made  here  to  specify  from  this  data.  For  the  following  dis¬ 
cussion,  is  not  a  critical  quantity  and  we  shall  assume,  as  in  the  theoretical 
curves  of  Figs.  14  to  18,  that  the  forward  lobe  of  the  phase  function  has  a  3  dB 
beamwidth  of  6°,  though  this  value  may  be  rather  small. 


(2)  The  extinction  cross  section  Oj  determines  the  attenuation  rate  at  short 
distances  into  the  random  medium  where  the  coherent  component  dominates.  With 
Equation  (40),  we  have  in  this  range: 

“T  -O-j-Z 


In  accordance  with  the  experimental  situation,  we  have  assumed  here  that 
°R  “  o  a  0.  The  relation  permits  one  to  determine  from  the  measured  range 
dependence  curves  of  Figs.  23  and  24.  Recalling  that  the  tree  spacing  is  *13  m, 
one  finds 

i  0.4  m  *  for  trees  in  leaf 
0.25  m 


oT  =  oA  •  o 


for  trees  without  leaves 


(43a) 


These  values  correspond  to  the  steep  initial  part  of  the  experimental  curves 
and  are  typical  for  both  nun-wave  frequencies  (20.8  and  07.6  GHz)  and  both  heights 
above  ground  (4  in  and  6  m).  However,  these  values  appear  to  be  very  high.  Since, 
in  a  forest,  all  scatter  objects  are  large  compared  to  a  mm  wavelength,  Oy  should 
be  close  to  the  high  frequency  limit,  which  is  obtained  forA  ->  0.  in  this  limit,  the 
extinction  cross  section  approaches  the  obstruction  cross  section  per  volume  of 
the  random  medium.  Thus,  according  to  Equation  (43a),  trees  in  leaf  would  obstruct 
*4  OX  of  the  planar  field  of  view  for  any  meter  of  vegetation  depth,  and  trees  with¬ 
out  leaves  would  lead  to  a  257!  obstruction  per  meter.  Although  values  of  this  mag¬ 
nitude  may  occur  locally,  it  is  very  unlikely  that  they  can  be  regarded  as  repre¬ 
sentative. 

A  possible  explanation  for  these  large  Oy  values  can  be  derived  from  the 
fact  that  in  a  forest  environment  all  scatter  objects  are  large  compared  to  a  mm 
wavelength.  Behind  large  obstructions  such  as  strong  branches,  deep  shadow  will 
exist*  and  since  the  measurements  were  pet  formed  over  a  path  through  maximum  vege¬ 
tation  density,  it  is  very  likely  that  the  receiving  antenna  was  operating  within 
such  shadow  zones.  The  measured  date  may  thus  indicate  a  higher  value  for  Oy  than 
one  would  obtain  in  the  average,  i.e,,  by  performing  similar  experiments  over  a 
number  of  parallel,  laterally  offset,  vegetation  paths  (including  the  maximum  den¬ 
sity  path)  and  by  averaging  the  powers  received  over  these  paths.  In  transport 
theory,  which  neglects  interference  effects,  this  averaging  process  is  implied. 

At.  large  vegetation  depths,  on  the  other  hand,  deep  shadow  zones  are  not  ex 
pected  to  exist  since  the  incoherent,  (muitiscattered)  component  is  well  developed 
in  this  region  and  will  "fill  in"  any  potential  shadow  zones.  Hence,  it  should 
bo  oossiblo  to  achieve  good  agreement  between  the  experimental  and  theoretical 
range  dependencies  at  large  r. 

Since  the  above  values  for  uy  appear  to  be  very  much  on  the  high  side, 
curve  fitting  of  the  theoretical  and  experimental  range  dependence  curves  was 
attempted  not  only  for  these  o,  values  but  for  moderately  and  strongly  reduced 
values  of  °y  as  well,  i.e.,  for 

0.3  m  ^  for  trees  in  leaf 

°T  =  '■  ,  (43b) 

0.2  m  for  trees  without  leaves 

c 

*  Note  that  we  are  concerned  here  with  the  range  of  small  vegetation  depths  where 
the  coherent  (highly  directional)  component  dominates  and  the  incoherent  com¬ 
ponent.  is  not  as  yea  significantly  developed. 


aT  --  S  (43b) 

~1 

0.05  m  for  trees  without  leaves 

i 

Using  these  values,  reasonable  agreement  should  be  achievable  at  moderate  and 
large  distances  into  woods.  At  short  distances,  discrepancies  are  unavoidable, 
since  the  assumed  theoretical  attenuation  rate  of  the  coherent  componer  differs 
from  the  measured  rate. 

(3)  The  albedo  W  and  the  ratio  a  of  the  forward  scattered  power  to  the 
total  scattered  power  were  adjusted  to  obtain  reasonable  overall  agreement  of  the 
calculated  and  measured  range  dependence  curves.  Agreement  was  established  simply 
by  trial  and  inspection.  Figs.  27  and  28  show  the  results  for  the  values  of 
Equations  (43a)  and  (43b),  respectively.  In  these  figures  the  received  power  in 
dB  (vegetation  loss)  is  plotted  as  a  function  of  the  normalized  distance  x  =  Ojl. 

The  solid  lines  show  the  theoretical  curves;  the  experimental  curves  (dashed  lines) 
are  those  of  Fig.  23  applying  to  a  propagation  path  at  a  height  of  4  m  above  ground. 

At  large  x,  good  overall  agreement  ‘  ^tween  the  theoretical  and  experimental 
range  dependence  is  obtained  in  particu  ,r  in  Fig.  28.  At  small  x,  the  agreement 
is  reasonable  for  trees  without  leaves;  for  trees  with  leaves,  the  measured  vege¬ 
tation  loss  is  significantly  higher  than  the  theoretical  loss,  in  particular  near 
the  sharp  bend  of  the  experimental  curve. 

Both  figures  indicate  that  the  parameters  W  and  a  must  be  chosen  close  to 
unity  to  achieve  reasonable  agreement,  i.e., 

W,  a  =-  0.9  to  0.55. 

For  the  small  o^.  values  of  Equation  (43b),  conclusive  results  were  difficult 
to  obtain  since  the  range  of  smail  x,  where  significant  deviations  between  the 
theoretical  and  experimental  curves  occur,  is  now  extended  and  the  range  of  large 
t,  where  good  agreement  can  be  achieved  in  principle,  is  correspondingly  reduced. 
However,  it  appears  that,  in  this  case,  the  conditions  on  W  and  a  are  somewhat  re¬ 
laxed  leading  to  values  of  0.8to0.95  for  W  and  0.6to0.9  for  a  (where  the  larger 
W  values  correspond  to  the  smaller  a  values  and  vice  versa). 

It  is  evident  that  these  results  are  preliminary  and  do  not  lead  to  definite 
conclusions  concerning  the  numerical  values  to  be  assigned  to  the  four  theoretical 
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FIGURE  27:  Comparison  of  Theoretical  and  Experimental  Range  Dependence 

The  experimental  curves  are  those  of  Fig.  23  and  apply  to  a  transmission 
path  at  4  m  height.  In  the  theoretical  curves,  the  parameters  W  and  a 
were  adjusted  for  good  agreement  with  the  experimental  curves.  The  ab- 
cissa  scale  refers  to  the  normalized  vegetation  depth  x  where  it  is 

assumed  that  oj  =  0.4  m-1  for  trees  in  leave  and  oj  =  0.25  m'l  for  trees 
without  leaves 


parameters.  It  appears,  however,  that  both  a  and  W  should  not  be  too  far  away  from 
unity  which  would  indicate  that  forward  scattering  is  a  significant  effect  for  mm- 
wave  propagation  in  woods  and  forests,  and  that  scattering  exceeds  absorption  by  a 
substantial  margin. 

V .  CONCLUSIONS 

A  theoretical  approach  for  the  study  of  mn-wave  propagation  in  woods  and  for¬ 
est  has  been  derived  and  compared  to  experimental  data  obtained  by  an  independent 
investigation  [1,  2).  The  theory  models  the  vegetation  environment  as  a  statisti¬ 
cally  homogeneous  random  medium  of  scatterers  characterized  by  a  scatter  function 
(phase  function)  with  a  narrow  forward  lobe  and  an  isotropic  background.  The  scalar 
transport  equation  has  been  used  to  determine  both  the  coherent  and  the  incoherent 
field  intensities  in  this  medium;  the  particular  problem  solved  is  that  of  a  vege¬ 
tation  halfspace  illuminated  by  a  plane  wave  Incident  from  the  air  halfspace.  The 
propagation  and  scatter  properties  of  the  vegetation  medium  are  characterized  by 
four  parameters  which,  in  principle,  can  be  determined  by  a  comparison  of  theoreti¬ 
cal  and  experimental  results.  Numerical  evaluations  have  provided  information  on 
the  range  dependence  of  mm-waves  in  woods  and  forests,  and  on  beam  broadening 
effects. 

Comparison  with  measured  data  has  shown  good  qualitative  agreement  between 
theory  and  experiments.  Both  indicate  that  the  attenuation  rate  of  mm  waves  in 
woods  and  forests  decreases  from  a  high  value  at  small  distances  into  the  medium 
to  a  lower  value  at  large  distances,  and  that  significant  beam  broadening  occurs, 
resulting  in  a  very  broad  beamwidth  at  large  vegetation  depths.  Reasonably  good 
quantitative  agreement  between  experimental  and  theoretical  results  can  be  estab¬ 
lished  at  large  vegetation  depths  (by  appropriate  adjustment  of  the  parameters  of 
the  theoretical  model).  But  the  experimental  data  indicates  a  higher  attenuation 
rate  at  short  distances  into  woods  than  can  be  accommodated  by  the  present  theory. 
The  comparison  between  theory  and  experiments  has  not  yielded  conclusive  quanti¬ 
tative  results  at  the  present  stage;  but  it  indicates  that  the  forest  environment 
at  mm  wavelength  is  a  strongly  scattering  medium  whose  scatter  cross  section  per 
volume  substantially  exceeds  the  absorption  cross  section  and  that  strong  forward 
scattering  occurs,  which  is  expected  since  all  scatter  objects  have  large  dimen¬ 
sions  compared  to  a  mm  wavelength. 

>r  study  is  needed  in  the  following  areas: 
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(1)  The  very  high  attenuation  rate  measured  at  short  distances  into  woods 
and  the  rather  sharp  transition  to  a  region  of  lower  attenuation  rate,  which 
occurs  after  transmission  through  a  few  trees,  requires  explanation. 

(2)  More  conclusive  (experimental  or  theoretical)  information  is  needed  on 
the  numerical  values  to  be  assigned  to  the  four  (averaged)  parameters  by  which  the 
theory  characterizes  the  forest  environment  at  mm  wavelength.  These  parameters  in¬ 
clude  the  scatter  cross  section  per  volume,  the  absorption  cross  section  per  volume, 
the  width  of  the  forward  lobe  of  the  scatter  characteristic  (phase  function),  and 
the  ratio  of  the  forward  scattered  power  to  the  total  scattered  power.  In  more 
general  terms,  the  characterization  of  the  forest  environment  at  mm  wavelengths  by 

a  set  of  suitably  defined,  measurable  electrical  parameters  is  a  very  important, 
but  up  to  now,  unresolved  problem.  The  theory  presented  here  is  a  first  attempt 
in  this  direction. 

(3)  The  present  theory  deals  with  the  simple  model  of  a  forest  halfspace 
with  a  planar  boundary  surface  illuminated  by  a  plane  wave.  A  more  realistic  model 
would  replace  the  forest  halfspace  by  a  two-layer  horizontal  slab  model  over  a 
ground  plane  with  suitably  defined  absorption  and  scattering  properties.  The  model 
is  indicated  in  Fig.  29.  Furthermore,  a  vector  transport  equation  formulation  of 
the  problem  is  needed  if  depolarization  effects  of  mm  waves  in  forests  are  to  be 
included  in  the  theory.  Analytical/numerical  difficulties  in  treating  a  model  of 
this  type,  however,  may  be  formidable. 
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APPENDIX  A 


Energy  Density  Associated  with  1^. 

The  integral 

2ir  i! 

U^Z)  -  f  I  I X(Z ;G  ,4.)  sin  0  d  0  d  4,  (A.l) 

0  b 

which  appears  on  the  right  side  of  Equation  (16b)  can  be  evaluated  in  closed  form. 
We  use  Equation  (16a)  and  perform  an  integration  in  G  ,4>  (over  ft  =  4ir)  on  both 
sides.  This  leads  to  the  following  differential  equation  for  (Z): 

°A  +  °S  , 

aU,  ■  Tos  r 1  (A. 2) 

cos  6p  +  (oA  +  (l-M)as)Uj  =  ao$  Sp  e  p 

We  have  made  use  here  of  the  relation 


2ir  ir 

—■r  y  /  q  (y)  Sin  0  d  0  d  4>  =  1 
o  0 

which  follows  from  Equations  (4)  and  (14a),  and  have  replaced  the  factor  cos  0  on 
the  right  side  of  Equation  (16a)  by  cos  0p.  The  approximation  is  discussed  in  the 
main  text. 


Equation  (A. 2),  a  linear  differential  equation  of  first  order,  can  be  solved 
in  closed  form.  With  the  boundary  condition  U.|(0)  =  0  which  takes  into  account 
that  Ij  (0;  0,<|>)  is  zero  in  the  range  of  forward  directions,  0  <  0  <  ^  (see  Equa¬ 
tion  (21a)),  and  negligible  in  the  range  of  backward  radiation,  ^  £  0  <  n,  we 
obtain: 


Uj(Z)  =  S  {  exp  (- 


°A  +  (1-“)05  7,  ,  °A  +  °S  7,  , 

cos  0  -)-exP  (■  cos  0 

p  p 


(A.  3) 


As  a  final  remark,  we  note  that  U^(Z)  is  proportional  to  the  energy  density 
associated  with  the  specific  intensity  Ij. 
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Solution  of  Recursion  Relation  for  Ij. 


The  coefficients  bn)  (O.'iO  of  the  specific  intensity  1^  are  determined  by  the 
recursion  relations  (Equations  (23)  of  main  text): 


bl  =  $  co tV  ‘  q(V 


2fi  IT 


K  (M>)  ■  t 


1  aW  1 


4ir  cos  0  m 

1- 


0  0 


m  =  2,3,. 


where  cosy  =  cos  ((!>-<(>")  sin  ft  sin  B'  +  cos  o  cos  o" 
cosy  =  cos  (d>-4>n)  sin  6  sin  e  +  cos  0  cos  0 

r  r  r  r 


For 


2 

q(Y)  =  (A~ )  e  yS 
\s 


’  AyS  <<:  ” 


these  equations  can  be  solved  by  assuming  that  bm  has  the  general  form 


»,„  «>•♦>  ■  B,»  ‘r/  e 

yin 


Y  2 

-<**-> 
ym 


with  Aym  <<  " 

m  =  1,2,— 
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(B.l) 


q(y)  sinO'dO-'dr  (B.2) 


(8.3) 


(B.  4) 


-  ;J 


V  \  \  •»  X 


■  Y'  %*  S'  S'  \*  s*  V 


Equation  (B.2)  requires  that 

/Yp  j  B 

0  2  «W  I  /  A  %2 

Bm  e  4*  cos  %  Wr J. 

Yin  '  yS  yin- 1 


'ii  2  y  2 

x  f  I  exp  {  -  (j*— )  -  (tt— P— )  )  s i n  O^dO'cty ' 
i  J  yS  ym-1 


(B.S) 


where  cos  yp  =  cos  (<f>p-<K)  sin  Op  sin  o  ' +  cos  0p  cos  0' 

In  evaluating  the  integral  on  the  right  side,  it  is  important  to  understand 
that  in  the  integrand  we  are  dealing  with  three  directions,  i.e.,  0,c)>  and  ^>4^ 
and  0',<p'\  For  the  purpose  of  the  integration,  the  directions  0,<|>  and  0p,<i>p  are 
fixed  while  Q'<p'  varies.  Note,  furthermore,  that  y  and  yp'  are  the  angles  which 
the  variable  direction  includes  with  the  fixed  directions  while  y  ,  the  angle  be¬ 
tween  0,<J>  and  0pt<f>p  is  a  constant;  see  Fig.  Bl. 

Since  it  is  assumed  that  the  "beamwidths"  A  -  and  A  are  small,  the  integral 

Y  j  )  ill 

will  have  appreciable  values  only  if  the  directions  0 ,q>  and  o  ,*  are  in  close 

r  r 

proximity  to  each  other  so  that  there  is  a  0",<K  range  in  which  both  y  and  yp'  are 
small  simultaneously.  Only  this  range  will  yield  a  significant  contribution  to  the 
integral.  To  determine  this  contribution,  we  introduce  a  local  coordinate  system 
on  the  unit  sphere  over  which  the  0  integration  is  performed.  This  coordinate 
system  is  indicated  in  Fig.  02  in  planar  projection.  The  coordinate  u  is  counted 
along  the  great  circle  through  0,f  and  0  ,<|i  ;  the  coordinate  v  is  counted  in  the 

r  r 

direction  normal  to  this  circle.  Both  u  and  v  are  measured  in  radians.  Evidently, 


r2  =  (■-£  +  u)2  +  v‘ 


Y"  a  (y£ 
Yp  1  2 


,2  J  Z 
u)  i  v 


and 


FIGURE  Bl:  Directions  0,  <j>  and  O',  and  O',  «(>'. 

For  purposes  of  performing  integration  in  Equation  (B.5). 
Direction  O',  ^ '  is  variable,  while  directions  0,  q»  and 
Op,  <|'p  are  fixed. 
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'  i  >  e.  y  £ 

/  /  exp  {  -  (-'-•)  -  (—2-)  }  sin  0'‘dO'd<j) 

oo  Ys  W-1 


-foo 

/exp < - +  (Y~~)2J  -  i4“z  +  r”2 1  v2  }dadv 
JL  \s  Vi  ^  A^ 


(B.6) 


•n 


A  2  ,  A  2 
Y$  +  Ym-1 


I  +  1 
A  2  +  A  2 
YS  Ym 


The  extension  of  the  range  of  integration  to  infinity  has  a  negligible  effect  on 
the  value  of  the  integral.  Equation  (B.5)  is  satisfied  provided 


1 


1 


1  ctW 


A  2  A  ,2,  A  .2 
•ym  yS  +  ym-1 

With  Equations  (B.l),  ( B . 3)  and  (B.4),  we  have 

$ 


D  g  _ _ 

m  ~  m-1  m  cos  0_ 


A  A 
Y1  “  YS 


R  ..  P  aW 

7/r'  COS  0. 


Hence,  wit!)  Equations  (B.6b) 


ym  V  ^  yS 


S  .  in 

B  -  (  aH  -  \ 

an  4tt  m!  cos  0 

P 


and  with  Equation  (B.4): 


,  ,  ,  . ,  5_p.  1  ,otW,m _ 4  9 

m  1  4 it  m/  p  ’  in- A  c 

P  Y> 


-TeL2 

”"V. 


(B .  7) 
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APPENDIX  C 


Asymptotic  Behavior  of  ^ 

The  asymptotic  behavior  of  for  large  t  is  determined  by  the  expansion  term 
associated  with  the  largest  eigenvalue,  S^,  which  is  the  only  eigenvalue  greater 
than  unity  and,  thus,  provides  the  slowest  attenuation  rate.  With  Equations  (25) 
and  (31) ,  we  have 


I2(x;p) 


—P.  e 
2ir  N 


SN  " 


*  ^ 

n=o  ,  _n 

1  "  s 

N 


for  x 


where  is  an  amplitude  coefficient  determined  by  the  boundary  condition  at  r  =  0 
and  is  the  largest  solution  of  the  characteristic  Equation  (32).  With  Equation 
(30b),  this  characteristic  equation  may  be  written  in  the  form 


M  E  _AW!.  =  1 

2  n-o  ,  ^n 

j  -  e 


(C.2) 


where  Apn  -  ^n+1  pn-l 

As  N  increases  the  sum  on  the  right  side  of  Equation  (C.l)  approaches  the  function 
(1-  r  )"*.  Since  Sw  >  1,  no  singularity  will  occur  and  we  have 


^(xlu)  -  J)  e  d  A  f  (y)  for  r  - 


WMCI  U  -  c  A 

f  M  -  i"_  a  A  -  (-  -N  1  *  Sa  = 

a'|J'  .  u  ’a  "  i  '  X  N  N  N  -  « 

s’  s 

*a  N 

Similarily,  the  sum  on  the  left  side  of  Equation  (C.2)  converges  towards  the  integ- 

Jl  , 

ral  /  (1-  -<jJ-  )  du  as  N  +  and  the  characteristic  equation  for  Sg  fakes  the 
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(C.4) 


The  quantities  S3,  A,  and  i\  are  shown  in  Figs.  7,  8  and  9  as  functions  of  W. 

a  a  a 

While  f  and  S,  are  fully  described  by  Equations  (C.3)  and  (C.4)  and  depend  on  W 

a  a  ^ 

only,  the  amplitude  Ag  depends  also  on  the  boundary  conditions  at  x  =  0  and  the 
direction  of  incidence  of  the  primary  plane  wave.  Fig.  9  implies  0p  =  0  and  bound¬ 
ary  condition  (21b),  i.e.,  ^((Kv)  =  0  for  0  <  li  <  1.  Numerical  evaluations  have 
shown  that  for  sufficiently  large  N  >  20,  the  highest  order  term  (k  -  N)  in  Equa¬ 
tion  (34)  approximates  Equation  (C.3)  closely. 

Equations  (C.3)  and  (C.4)  can  also  be  obtained  directly  by  rigorous  solution 
of  Equation  (19b)  with  the  forcing  term  suppressed.  This  term  decreases  with 
exp  (-x/pp)  and  is  insignificant  at  large  c.  This  approach,  however,  leaves  the 
asymptotic  amplitude  Aa  undetermined. 

a 


